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Abstract. The Be star phenomenon is related to fast rotation, although the cause of this fast 
rotation is not yet clearly established. The basic effects of fast rotation on the stellar structure 
are reviewed: oblateness, mixing, anisotropic winds. The processes governing the evolution of 
the equatorial velocity of a single star (transport mechanisms and mass loss) are presented, as 
well as their metallicity dependence. The theoretical results are compared to observations of B 
and Be stars in the Galaxy and the Magellanic Clouds. 
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1. Introduction 

Why care about rotation? Just becaus e stars do rotate! A look at the velocity distribu- 
tion established bv lHuang fc Giej (2006*) reveals that the peak in the probability density 
occurs at Vcq — 200 km/s, which represents a ratio v/vait — 0.5 — 0.6, i.e. a substantial 
fraction of the keplerian velocity. 

The effects of rotation in stars were studied since the works of von Zeipel and Eddington 
in the years 1924-1 925. In the end of the 60s, they were included in polytropic or sinipli 
fied s t ellar models (IRoxburgh et al.l 19651: [Roxburgh fc Strittmattenl 19661 : iFaulkner et al 



1968t iKippenhahn fc Thomaslll970l : lEndal fc Sofialll976[ l. AbouTsO years later, stellar 
mode ls became more sophisticated and also benefited from the inclusion of rotational ef- 
fects (IPinsonneault et aL 19891 : Deupred 1990l : Fliegner fc Langer 1994; Chabover et al 
1995t lMevnetlll996l) . 



Since th e end of the 90s, niore or less extended grids of rotating stellar models were 



comp uted (jLanger et al.lll997t iMevnet fc Maedeijll997t ISiess fc Liviolll997t iHeger et al 



l200C/). Those grids showed that the inclusion of the effects of rotation improved the 

adequation between models and massive stars observations in many aspects: 

» t he surface abundances of light elements (jHeger fc Langer l200nl: iMevnet fc Maedeii 
200n[) : 



20001) : 



the predicted surface velocities in clusters ( Martavan et al.|[2006al lbl: Mevnet fc Maedei 



the blue- to red-supergiants number ratio in the SMC (IMa eder fc Mevnet '200% 



• t he WR populations number with metallicity (jMevnet fc Maeder..2003.,2005; ,Vink fc de Koter 
200i) 



the rotation rates of pulsars (when strong coupling is assumed. Heger e t al.''2005): 



20091) 



the SN types and GRB progenitors (Mevnet fc Maederii2005i : iYoon et al.l2006i;,Georgy et al 
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Roughly summarised, rotating stars are expected to present 
(a) a modified gravity: 

• the surface characteristics become dependent on the colatitude considered; 

• there is a mass loss enhancement and anisotropy; 

(&) chemical species and angular momentum transport mechanisms: 

• the behaviour on the HR diagram is modified; 

• the nucleosynthesis is altered; 

• there is a surface enrichment; 

• the mass loss is modified; 

• the rotation profile evolves during the star life, becoming steeper (when no mag- 
netic fields are considered). 

We will go through the different points in this review. 



2. Surface characteristics 

2.1. Gravity and shape 

Because of rotation, the effective gravity is modified, becoming a function of the rotation 
velocity and of the colatitude 9 of the star: 

5off = 9cs{^, 0) = (^-^ + sin^ 9^ eV 
4- ri^r sin^ cos 6* eg. 

We immediately see that at the pole {0 = 0°) the effective gravity is just the gravitation 
acceleration —GMjr^. At the equator (9 = 90°) the centrifugal force adds a sustaining 
term sin^ 9. In these two cases, the effective gravity is still radial, while at interme- 
diate 6*, the term sin0 co%9 does not vanish and implies that the effective gravity is 
no more radial. 

In the frame of the Roche model, the maximal oblateness allowed when the star rotates 
at the critical velocitjU is i?oq,crit — 1-5 i?poi,crit. 

Recently, interferometry has allowed to determine the de formation of fast rotating 
stars. A first evaluation of the oblateness of Achernar by iDomiciano de Souza et al] 
(l2003h showed a larger ratio i?oq/-Rpoi than th e one allowed in the Roche model. How- 



ever, more recent observations of the same star (|Vinicius et al.ll20Q6t ICarciofi et al.ll2008l ) 
have revised this ratio and found a lower value, more compatible with these theoretical 
expectations. 

2.2. Flux and effective temperature 



Since the flux is re lated to the effective gravity ()von Zeipell Il924t lOwocki et al.l 11996 



1998HMaedeJll999l) . it becomes also dependent on colatitude 



with M* = Al — rt-nGp ) ^'^^ so-called reduced mass which takes into account the 
reduction of the gravitational potential by rotation (pm is the mean density inside the 
considered isobar). 

The Stefan-Boltzmann law F = aT'^ implies a dependence on the colatitude for the 

t The critical velocity is reached when g^B vanishes because the centrifugal force at the equator 
counterbalances the gravity exactly (see Section [3. 2 |l . 
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Figure 1. Variation of tiie effective temperat ure with the colatitu de for various initial rotation 

rates (figure from lEkstrom etlll I2008D . 



effective temperature as well: 



iiraGM 



7 ges{i^, t 



1/4 



or 



Recent interferom etric observations, like for example iMonnier et al. ] (|2007t) on Altair 
Zhao et al. I (l20Q9l) on a Cephei and a Ophiuchi, have provided a possibility to test 
this relation. Let's have a close look at Altair: the rotation rate fi/ficrit of this star is 
evaluated betvifeen 0.90 and 0.92, and the temperature difference betvifeen the pole and 
the equator is betvifeen 1.19 and 1.32. For such a rotation rate, the theoretical models (cf. 
Fig. [T]) predict a difference of 1.26-1.33, in good agreement with the observational value. 



2.3. Mass loss 

Accor ding to the works of lOwocki fc GavlevI |l997l ). lMaeder fc Mevnetl (|2000l) and lPetrenz fc Puld 
(|2000l) . rotation enhances the mass loss by a factor: 



M{n) 

M(0) 



(1 - TEdd) 



1 - 



27T G p„ 



Tsdd 



where Fsdd is the Eddington factor, i.e. the ratio of the luminosity of the star to the 



Eddington luminosity iEdd = 



with Kg the electron-scattering opacity. 



4ir c G A/ ' - . 11 _ I . - - 

It also changes the geometry of the mass flux (jMaeddl2Q0al2Q09f), which is no longer 



constant on the whole surface of the star. The mass loss by surface unit at a given latitude 
9 follows the relation: 



dM(6>) 
do- 
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Figure 2. Left: evolution of the equatorial velocity (normalised to the initial velocity) during 
the main sequence for different mass domains. Right: evolution of the equatorial velocity, the 
critical velocity and the ratio v/vait during the main sequence. 



where A{a,k) — (fca)^/" i^'TT')^^ ^ function of the force multiplier parameters 

which characterise the stellar opacity. 

If we consider that A = est, a star rotating for example at il/ricrit = 0.95 will have 
Af(pol) — 3.2M(eq). However, if one accounts for the change of the force multiplier 
parameters, i.e. a variation in the opacity regime caused by the drop of Tcff at the equator, 
the equatorial mass loss can be enhanced at a given point, driving the formation of a 
decretion disc (IOwockill2004() . 

On the observations side, interferometry again sheds a new light in this topic, showing 
featu r es that can b e interpreted as polar enhanced winds (* Kervella fc Domiciano de Souzal 
20061 iMeilland et al . 2007), as well as discs around active stars (^Meilland et al. 2 007; 
Schaefer et al.ll2010l) . 

The geometry of the wind may leave an imprint on the circumstellar medium. There 
are some indications of asymmetry detecte d in spectropolarime try observations of some 
supernovae (see for example SN 2007rt bv lTrundle et al.ll2009t ). This aspect is actually 
under study with 2- and 3D simulations and will be the subject of a future paper (Walder 
et al. in prep). 



3. Rotational evolution 

3.1. Two competing processes 

The evolution of the surface velocity of a star is the result of the competition between 
two processes: 

(a) the mass loss, which removes angular momentum at the surface and thus decel- 
erates the rotational velocity; 

(6) the transport of angular momentum inside the star, which brings some internal 
angular momentum to the surface (mainly through the meridional circulation) and may 
counterbalance the loss by the winds. 

Both processes are dependent on the mass of the star: the more massive the star, the 
stronger winds and larger meridional currents it experiences (see Fig. [21 left). At solar 
metallicity typically, around 20 Mq and above, the winds contribution wins so the star 
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Figure 3. Evolution of tlie v/vc-cn ratio during tlie main sequence for 20 Mq models at 

Z = 0.020 and Z = 0.002. 



decelerates during the main sequence (MS). Around 9 Mq, both processes are in equi- 
librium so the surface velocity remains almost constant during the largest part of the 
MS. Around 3 Mq, there are almost no mass loss through winds, but the core-envelope 
coupling through the meridional circulation is s o weak that the surf ace velocity decreases 
during the MS. Note that recent observations bv lHuang et al. ( 2010l see also Poster Sl-05 
in these proceedings) show a quicker deceleration of the less massive of the B stars than 
predicted by current theoretical models, more compatible with a regime of local angular 
momentum conservation Q.E? = est. 



3.2. Critical limit 

One thing is the evolution of the surface velocity, another thing is the evolution of the 

ratio to the critical limit Vcrit — During the MS, the mass may be reduced 

by mass loss mechanisms, and the radius steadily inflates, so the critical limit drops, as 
shown in Fig. [2l right. The result is that although the surface may decelerates, the star 
may encounter the critical limit at a moment in its main sequence lifetime. Comparing 
the left panels of Fig. |3] and 21 we see that the conditions to reach the critical limit during 
the MS are met if the star has not a too high mass, and also if it is not a too slow rotator 
at birth. For each mass domain, there is a minimal initial u/ucrit ratio allowing for the 
reaching of the critical limit. 

Once at the critical limit, the star remains close to it. It experiences phases of mass 
ejection (slowing the surface below the critical value) followed by quiescent phases, during 
which it slowly re-accelerates toward the critical limit. 

It is possible now for theoretical models to evaluate the mass ejected in the form of a 
disc (see Poster Sl-03 by Georgy et al. and Sl-06 by Krticka et al. in these proceedings). 
This 'mechanical' mass lo ss seems to occur at a lower ra te than the one that can be 
measured around Be stars ( Rinehart et al. 1999t Stee 2003 ). However, in the models, the 



mechanical mass loss is averaged on a much longer timestep than the period of mass 
ejection observed in Be stars. The instantaneous mass-loss rate is expected to be higher 
than the average one given by the models. 
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Figure 4. Evolution of the v/vcrit ratio during the main sequence for 9 Mq models at 

Z = 0.020 and Z = 0.002. 



3.3. Metallicity effects 

Both the radiative mass loss rate and the strength of the meridional circulation depend on 
the metallicity. A low metallicity reduces both mechanisms: there are less metal lines to 
interact with the photons and form a radiatively driven wind and the greater compactness 
of the star reduces the amplitude of the meridional circulation. 

The net effect on the evolution of the equatorial velocity depends strongly on the mass 
domain considered: for the most massive stars, the reduction of the winds is the dominant 
effect, so the stars reach more easily the critical limit (see Fig. On the contrary, for 
less massive stars for which the winds are not too strong anyway, the reduction of the 
meridional circulation is dominant, so the stars have more difficulties to reach the critical 
hmit (see Fig. HJ. 



3.4. Be phenomenon 

Fast rotation is supposed to be linked with the Be phenomenon. It is however not yet 
clear whether it may explain by itself the origin of the eq uatorial disc observed around 
Be stars. Obs ervations show a ratio u/ucrit ~ 0.70 — 0.80 (jPorteiillQQS IChauville et al 



2001 



Tvcner et al. 20051 ). Ma ybe this ratio is t rue, b ut there may also be alternative 

explanations: first, as shown bv lTownsend et aL I (l2004l) . there is a saturation effect in the 
widening of the lines by rotation, so the true velocity may be underestimated. Second, the 
Be phenomenon seems to consist of ejection phases followed by quiescent phases during 
which the stars may rotate with a lower rate. 

Observati ons show a me tallicity tren d in the appearance of the Be phenomenon (see 
for example Maeder et al. |1999 or Wisniewski fc Biorkman 20061 ). Two scenarios are 
evoked to explain the f ast rotation of the Be stars: the binary channel (see for example 
McSwain fc Gies 20051 ) (where the fast rotation arises from the accret ion of angular 



mom entum from the mass-donor companion), or the single star channel ([Ekstrom et al 



20081 ) (where the surface acceleration is a natural evolution due to the core-envelope 



coupling). Any such scenario sho uld present this metallicity trend to be valid. To study 
the single star evolution scenario, Ekstrom et al.l (.2008) computed 112 stellar models (4 
masses, 4 metallicities and 7 rotation rates) with the Geneva code. To get populations 
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n umbers, we conyolved the models with a Salpeter IMF, using the velocity distribution 
of lHuang fc Gie^ (l2006l ). 

A good agreement between the theoretical population ratio Be/(B+Be) and the ob- 
served one is obtained modulo two adjustments: 

(a) the Be phenomenon should appear already at 70% of the critical velocity; 

(&) the velocity distribution should count more fast rotators at birth in low-metallicity 
environments. 

Point (a) can be sustained by several mechanisms. For example, sub-surface convective 
motions have been shown to be able to give the n eeded impulse to eject matter from the 
surface ([Maeder et al. 2008 : Cantiello et al.ll2009l) . as also would non-radial pulsations do 
( McSwain et al.ll20Q8t lOwockilbOO^ . The binarity of the star could play a role, here not 



by spinnin g-up the star but by adding a gravitational pull allowing to launch the matter 
in the disc ( Kervella et al1l2008l) . Point (b) presents contradictory observational supports. 
While Martavan et al.l ( 2007) do find a higher me an fl /flcrit ratio in the SMC compared 
to the LMC and the Galaxv. lPmnv fc GiesI (|2009( ) draw opposite conclusions, finding no 
clear evidence for any difference in the velocity distribution at different metallicities. 

Note that the stellar models have been computed without the effects of magnetic fields, 
the inclusion of which could change the picture. It may be that taking into account the 
strong core-envelope coupling brought by the magnetic fields relieves the theoretical 
prediction from the two aforementioned adjustments. 

3.5. LBVs 

Until now, the critical velocity we have considered was determined only by the centrifugal 
contribution (/lot against .ggrav (fi-limit). However, a critical point appears whenever 
9 tot = ffgrav + 5rot + 5rad = 0. If grot = Or is negligible, the star may meet the classical 
Eddington- limit when ^rad = 5grav If .9rad is negligible, the star may encounter the ri-limit 
when grot = ggrav- this describes the first critical velocity considered until now: 



Wcrit,l 



12 GM 

3 -Rpol,crit 



with i?poi,crit the polar radius when the star is at the critical limit. 

If all three terms are significant, the star may meet the riF-limit when ^tot 
means that for a given fi, there is a maximum luminosity given by: 



0. It 



AttcG M , ^ 



2tT G pn 



with K the total opacity. Inversely, depending on the Eddington factor of the star, a 
second critical velocity can be defined: 



i-r„ 



Wcrit,2 



7?2 

"'pol,crit 



where the quantity V'{io) = 



"'^'-^ V'{u;) 

is the ratio of the actual volume of a star with 



rotation uj — Q/Qcnt to the volume of a sphere of radius -Rpoi,crit. 

Note that for Fsdd < 0.639, z;crit,2 is not defined. Above this value, Wcrit,2 becomes 
smaller than Ucrit.i, so the star enc ounters the r^F-limit be fore the ri-limit. This could 
be the case of some known LBVs (jCroh et al.l l2006l l2009t ). that present both a high 
Eddington factor and a high rotation rate. With a Fsdd — 0.8, they are probably meeting 
their second critical velocity already with J7/f2ciit — 0.85 — 0.9. 
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Figure 5. Left: evolution of the surface N/C ratio (normalised to the initial one) during the 
main sequence (expressed as the central mass fraction of hydrogen Xc) for 9 Mq models at 
Z = 0.020 with varying the initial ratio rj/ficrit- Right: same as left, but for different mass 
domains. 

4. Internal mixing 

4.1. Diffusion of chemical species 

The radiative zones inside a star are supposed to rotate differentially. This causes shear 
turbulence which in turn drives some turbulent mixing of chemical species. The higher 
the rotation rate, the stronger is the mixing (see Fig. \E\ left). 

We can evaluate the diffusion time: Tmix — %■ with D the diffusion coefhcient. With a 
relation between the radius and the mass, we get Tmix cx M~^'^ while the main sequence 
lifetime is tms oc A/^'^-^. We see that the more massive the star, the stronger is the 
mixing, leading to a surface enrichment already on the main sequence (see Fig.[5l right). 

4.2. Metallicity effects 

Low-metallicity stars are more compact than higher-metallicity ones. The meridional 
currents are less efficient, so the i7-profile inside the star is steeper. The mixing is thus 
stronger, while the mixing time is shorter, so we expect a strong surface enrichment on 
the main sequence. 

4.3. Surface abundances 

The surface abundances are expected to be modified by mixing. A good signature of 
mixing is the C, N, and O abundances: their respective ratios are expected to change, 
but their sum C+N+0 is supposed to remain constant. Only at very low Z, primary 
nitrogen production ma y occur, leading to a n et metallicity increase. A careful abun- 
dance analysis made by Przvbilla et al. ( 201Cl() in the solar neighbourhood shows that 
the observed mixing follows very well the trend expected from CNO nuclear reactions 
(see their contribution in the same proceedings). The adequation with stellar models is 
also good, but the observations are not very constraining yet. 

It is very important to keep in mind that the mixing is a function of the rotation 
rate, of course, but also (as shown above) of the mass, the metallicity and the age, 
as well as other characteristics as the binarity for example. The study of the rotational 
mixing in stars should imply a separate analysis, in order to discriminate between several 
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parameters. For example, while iHunter et alj (|2Q08l ). studying field and clusters stars in 
th e SMC, find many ou tliers from the expected path in the N/C vs Urot, a re-analysis 
by iMaeder et al. ( 20091 ) of a sub-group of the same stars, restrained to a narrow mass 
range and belonging to the same cluster (i.e. having the same age), shows that most of 
the outliers shifted back to the expected path. Most of the outliers left cannot be used 
to test the rotation-induced mixing, because they are evolved stars or known binaries. 
For the very few remaining, we indeed need another explanation than pure rotational 
mixing. There are several possibilities that shall be the subject of further studies, as the 
influence of magnetic fields inside the star or magnetic braking at the surface. However, 
the need for a yet unknown process is evoked (see the contribution of Ines Brott in the 
same proceedings). 



5. Conclusion 

Rotation is linked to many very interesting types of objects as Be stars or LBVs. 
The inclusion of its effects in theoretical stellar models is essential to understand those 
objects and their direct environment. While rotating models are in a better adequation 
with observations, there is still a lot of work to do in order to improve our understanding 
of the physical processes at work. 

Theoretical and numerical developments are continuously on-going, but theory alone 
is like a car without a driver. Observations are highly needed: 

• larger surveys of well identified objects would allow to make a separate analysis of 
sub-groups in masses, ages, with a significant statistics; 

• a larger number of objects studied by interferometric measurements would allow to 
put constraints on the surface characteristics of the models; 

• the progress in asteroseismology should provide constraints on the internal structure 
of the stars; 

• precise observations of the circumstellar environment could help to constrain the 
mass loss mechanisms, which are a key ingredient in massive stars. 
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Discussion 

S. Owocki : I would like to emphasize an important difference between reaching the 
two critical rotation speeds identified in the Maeder & Meynet analysis. The first critical 
speed applies to low luminosity stars like Be stars, and should lead to a circumstellar 
decretion disk that is ejected mechanically at the equator. But for the second critical 
speed, which is modified by the radiative forces associated with T > 0.65, the mass loss 
should be mainly over the poles, not equator. In effect, the rapid rotation and associated 
equatorial gravity darkening forces the stellar luminosity to emerge over a smaller surface 
area over the poles, so that even if F < 1 the local flux over the poles can exceed the 
Eddington value, leading to a radiatively driven, bipolar, prolate mass loss. In short, in 
low luminosity, F << 1 Be stars, critical rotation leads to the observed equatorial disk, 
while in F < 1 LBVs it leads to bipolar nebulae. 

O. Chesneau: a comment to Stan: If the mass loss of LBVs is essentially prolate, and 
directed toward the pole, how can you explain why the environment of LBVs is dominated 
by ring-like structures (see comment of K. Weis)? 

S. OwocKi: Well, in slowly rotating LBVs you should get spherical mass loss that will 
likely appear ring-like. Perhaps Kerstin can comment on how relatively common these 
rings are vs. bipolar LBVs. 

K. Weis: As Stan already mentioned, LBVs create polar winds which lead to larger 
number of bipolar nebulae and not only spherical ring nebulae. See my talk on Thursday. 



